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A Kinematic Model of 
the Human Spine and Torso 



Gary Monheit and Norman L Badler 
University of Pennsylvania 



Improving human mo- 
tion models requites a 
more realistic model of 
the body, including a 
flexible torso and spine. 
Defining spine attri- 
butes and input parame- 
ters helps achieve 
inverse kinematic con- 
trol of the spine. 
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searchers have siudicd human figure models in computer graphics 
almost since the introduction or the medium. "Jit rough the Inst dozen years 
or so. the structure, flexibility, and fidelity of human models have increased 
dramatically— from the wireframe stick figure, through simple polyhedral 
models, to curved surfaces, and even finite-clement models. Computer 
graphics modelers have tried to maximize detail and realism while maintain- 
ing a reasonable overall display cost. 

The same issue pertains to control: Improving motion realism requires 
many decrees of freedom in the body linkage, and such redundancy strains 
effective and intuitively useful control methods. We can cither simplify 
control by simplifying the model, thereby risking unrealistic movements, or 
complicate control with a complex model and hope the resulting motions 
appear more natural. The recent history of computer animation of human 
figures has focused on (he quest to move the technology front the former 
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Figure 1. Contoured body bent ton target slate vector: (flexion, 
axful, Infernl) e (-52 degrees, $H degrees ~7 degrees). 



situation towards the. latter while simultaneously forcing the. 
eumrol complexity into algorithms rather than skilled manual 
manipulation. 

'Phis point of view motivates our efforts in human figure mod- 
cling and aninwtion. a?> well as ihosc of several other groups. In 
particular, notable algorithms for grenier animation power 
have addressed kinematics, dynamics, inverse kinematics, 
available torque, global optimization, locomotion, deforma- 
tion, and RCKtural and directional control. 1 Throughout this 
ranpe of studies, however, the human models themselves 
tended to be raibei simplified versionsof real human flexibility. 
In the early 1 l As0s we warned that increased realism in the 
models would demand ever more accurate and complicated 
motion control. Now that the control regimes are improving, it 
is time to re turn to the human models and ask if we must 
reevaluate their structures to take advantage of algorithmic 
improvements. When we considered this question, we deter- 
mined that a more accurate model of the human spine and torso 
would be essential to further realism in human motion. 

Although many models have appeared to have a flexible 
lovso* they have been computer constructions of the surface 
shape manipulated by skilled animators. 2 YVc needed a torso 
suitable for animation that also satisfied our research project 
requirements for anthropometric seal ability/ Thus, wc could 
not accept n single model of fixed proportions when we would 
potentially need to model vast differences among human body 
types. 

Snakes exhibit a similar type Of flexible figure, 4 but ihe an- 
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Figure 2, Contoured body ben I lo a target state vector* (flexion, 
axial, lateral) c (-] degree, -60 degrees, -6 degrees)* 




Figure 3. Contoured body bent to a Hii-rcI stale vector: (flexion, 
axial, Jalcrat) = (3 decrees, -42 degrees, 23 degrees). 



thropometry issues do not arise. Miller based his animation 
approach on dynamics. Because humans do not need lo 
loeomote with their torsos, we deemed a kinematics model 

f fit 7 . IS Computer Graphics & Applications 
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Hgnre 4, Contoured horly bent to a target stale veiior; (fluxion, 
axial, lateral) = (-10 degrees, 10R degrees, -45 degrees). 



Figure 5. Contoured body bent to a target state vector (flexion, 
axial, lateral) * (J8 degrees, 108 degrees, -68 decrees). 




Figure*; CoibfoOfc'd body Wit loa tafgetstate viclfjfa (flexion, 
axlwl, lateral) « degr ecs, 39 degrees, 6 degrees), 



adequate. On the other hand, Zeltterand Stredney's "George" 
skeleton model had n derailed vertebral column, bur it was not 
articulated hot did it bend during Kinematic animation. 3 
WiUmen simulated limited neck vertebral motion in the 
wi^itat planed 

The default polyhedral figures used in "Jack*' (the software 
system developed at the University of Pennsylvania for human 

March 1991 



figure modeling and manipulation 7 ) lacked much detail it) the 
h uman torsos. The graphically displayed bodies appeared to be 
impersonating robots with stiff backs, bending only from the 
waist. A five-segment torso with more possibliiies of articula- 
tion was created a year ago, But the back was modeled without 
any curves, contrary to actual human anatomy, if the spine 
were realistically modeled, then the torso, a vessel connected 
and totally dependent on the spine, could be viewed and manip- 
ulated interactively. So one of us (Monheh) undertook the 
development of a far more satisfactory and highly flexible ver- 
, tebral model of the spine and its associated torso shape. 

The conceptual model of the spinal column is derived from 
medical data and heuristics related to human kinesiology. The 
spine Is a collection of vertebrae connected by ligaments, small 
muscles, vertebral joints (called processes), and intervertebral 
disks. 8 Nature has designed the spine for* 10 

1. Support of the body Is weight 

2. Stability of (he torso 

3. Flexibility of motion 

4. Protection of the spinnleord 

'lite spine moves as a scries of vertebrae connected by depen- 
demijohns," meaning that it is impossible to isolate movement 
of one vertebral joint from the surrounding vertebrae. Muscle 
groups of the head, neck, abdomen, nnd back initiate move- 
ment of the spine, and the interconnecting ligaments allow tho 
movement of neighboring vertebrae. 8,12 
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Figures 7-12. Wireframe polygonal body with viable spine 
moves through a rolling sequence. Zero Interpolation flag Is set 



The following sections describe both Ihe model and the rep- 
resentation of ibe model within Jack. 



Spinal motion 

The model of spinal nioiion is bused on 

• Anatomy of the physical vertebrae and disks 

• Range ol movement of each vertebra 

• I -If feci of llic surrounding ligaments and muscles 

Anatomy of vertebrae and disks 

The spina I column consists ofttvcrlc brae organized in to five 
regions 8 : 

• Cervical 

• ITiorncic 

• lAimbur 

• Sacral 

• Coccyx 

The vertebrae are labeled by medical convention: CI -C7,T1 
• T12, 1-1 - K5, and SI - S5. The regions listed above axe in 
vertical descending order. Which regions should wc. consider 
p;irl of the torso? The cervical spine lies within the neck. Also, 
I he sacrum and coccyx con tufn vertebrae fixed through fusion. 1 ' 
So, the mobile part of the torso includes the 12 thoracic arid 5 




lumbar vertebrae. Therefore, the torso model includes 17 ver- 
tebrae ond IS joints of movement. 

Each vertebra is uuiquely sized and shaped, but all vertebrae 
contain a columnar body and an arch. The body is relatively 
large and cyllndrlcaLsupportlng most of the weight of the en tire 
spine. The vertebral bodies increase gradually in size from the 
cervical to the lumbar region/'' 

The arch supports seven processes: four articular, two trans- 
verse, and one spinous," 'llic processes are bony protrusions on 

IEEE Computer Graphics & Applications 
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Intervertebral disks separate each vertebra. A disk lifts three 



parts": 

• nuclei!* pulposiis; the sphere in the centen consisting of HS 
percent water 

• annuhjs flbrosus: the fibers running as concentric cylinders 
around the nucleus 

• cartilaginous plates: a tbin wall separating the disk from 
the vertebral body 

'the disk changes shape as l he neighboring vertebrae bend. 
But, since the nucleus Is 85 percent water, very little compres- 
sion occurs. The disk can bulge out. spherically as force Is ap- 
plied to the columnar body above or below. Therefore, overall 
the disk does not function as a spring, but as a defonuable 
cylindrical separation between vertebrae. This supports the 
theory that the vertebrae do not slide , but rotate around an 



the vertebra that aid and limit the vertebral motion. The trans- 
verse and spinous processes serve as levers for both muscles and 
ligaments.* The articular processes provide a jVuit facet for ihe 
joint between successive vertebral arches. These processes, due 
to their geometry, cause the vertebrae to rotate with three de- 
grees of freedom. Ligaments and small muscles span successive 
vertebral processes, giving the spinal column its stability. Be- 
cause of this strong interco^incctivity, we model spinal move- 
ment as Interdependent movement of neighboring joints. 

March J 99 1 



axis. 



Vertebral range of movement 

Vertebral movement is limited by the relative size of the 
disks, the attached ligaments, and the shape and slant of the 
processes and facet joints. Statistics for joint limits between 
each sucessive vertebra have been recorded and compiled," 
Also, the spine has a natural shape at rest position. We use the 
initial joint position of each vertebra as input to the model. 

The range of movement differs for each region of the spine. 
For instance, the optimum movement of the lumbar region is 
flexion or extension. 7^e thoracic area easily moves laterally, 
whtlemthesagittaIpla7ieflexionandextensionareIimiied,Thc 
cervical area is very flexible for both axial twisting and lateral 
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Table L Rotation of the spine. j 


Flexion/extension 


Forward and 
backward bending 


Rotation around 
the* axis 


Axial rotation 


Twisting 


Rotation around 
the y axis 


Lateral bending 


Side bending 


Rotation around 
the z axis 



bending. Hie joint limits for each region affect how much that 
joint can participate in any given movement. The posture of the 
torso results from the specialization of the spinal regions. 11 

Effect of ligaments and muscles 

The vertebrae are interconnected by a complex web of sur- 
rounding ligaments and muscles. If the force initiated by a mus- 
cle group is applied at one joint, that joint moves, as do the 
neighboring Joints to a lesser degree. Some joints farther away 
might not be affected by the initiator joint's movement. 

It is possible, by simultaneous contractions of extensor and 
flexor muscles around the spinal column, to deactivate joints 
not initiating the movement. 13 Depending on the force of these 
resisting muscles, the joints on or near the joint closest to the 
resistor will move less than they would if the resisting force had 
not been applied The final position of the spine is a function of 
the initiator force, the resisting muscle, and the resistance. 

Input parameters 

We modeled the spine as a black box with an initial state, 
input parameters, and an output state. To initiate movement of 
the spine, we introduced several input parameters: 

♦ Joint range "from" and "to": Within the total number of 
joints in the spine, any nonempty contiguous subset of vertebral 
joints can be specified by two joint indices. These joints indicate 
which part of the spine is active in movement For example, the 
user specifies movement in the range between T5 and T10. All 
other joints are frozen during the movement 

• initiator joint: The joint where movement begins, usually 
the joint with greatest motion, 

• resistor Joint The joint that resists the movement This can 
bo equated to a muscle that contracts and tries to keep part of 
the spine immobile. 

♦ resistance: The amount of resistance provided by the resis- 
torjoint. 

* spine target position: A 3D vector describing the target 
position after rotation around the x, y, and z axes. The target 
position is the sum of all joint position vectors in the spine after 
movement succeeds. 

♦ te.ro interpolation: A M yes" indicates that movement is 
interpolated through the Joint rest position. A "no H indicates 
that only the joint limits are used to interpolate movement. 

34 



Vertebral joint movement 

The joint between each vertebra has three degrees of rota- 
tion, The spine will move toward the target position by rotating 
around the three possible axes, 11 as shown in Table 1. 

Hie position of the flexion rotational axis for each vertebral 
joint has been measured from cadavers. It is not equidistant to 
the two adjacent vertebrae, but lies closer to the bottom verte- 
bra. 11 The origin of the axis of movement determines how the 
vertebrae move. When the torso is modeled on the spine, the 
axis also directly determines how the torso changes shape. 

Elongation and compression are absent from the model. The 
hydrophilic intervertebral disk, when submitted to prolonged 
compression, induces a slight decrease in height due to fluid 
leakage. Conversely, after a lon$ period of rest or zero gravity, 
the spine elongates by maximum filling of the nucleus puiposus 
(at the center of the disk). 11 Dehydration during a day's activity 
can result in a loss of height of 2 centimeters in an adult In any 
short duration of movement the disk is essentially incompress- 
ible, and therefore elongation is imperceptible. 14 

Shearing or sliding (translational movements) of the verte- 
brae would lead to variation in the intervertebral separation, 
The mechanics of the intervertebral disk would not allow this. 1 1 
Therefore, we assume that for norma] activities the three de- 
grees of rotational movement are the only ones possible for 
each vertebral joint. 

Spine database 

Any human figure can have a wide variety of torso shapes. 
Also, each person has a different degree of flexibility and range 
of movement. To model the position and shape changes of an 
individual's spine, we designed a database for creating a unique 
set of features for the spine and torso. Medical data provides the 
database elements for an average person. 11 

Database elements include 

• Size of vertebra— x fi y,z dimension 

« Intervertebral disk size— separation between vertebrae 

• Joint limits— three rotations, two limits per rotation 

• Joint rest position — initial joint position of the spine 

Applying the motion model 

Consider a stationary spine and its attributes. Each vertebra 
has a current position defined by the joint position between 
each vertebra for each of the three degrees of rotation. The 
spinal database also includes definitions of joint rest positions 
and six joint limits for every joint If each attribute is summed 
up for all joints, then 3D vectors are defined for current posi- 
tion, joint rest position, and two joint limits for the global spine. 
The target position^-the 3D vector sum of final joint posi- 
tions — is supplied as an input parameter. Movement towards 
the target position is either bending or unbending, meaning 
either towards the joint limits or towards the spine's rest posi- 
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tion , Motion is defined as an interpolation between the current 
position a nd either the spine's position of maximum limit or the 
splne'srest position. 

Three rotations are calculated independently and then 
merged into one. For example, a 3D orientation vector (such as 
bend forward 45 degrees, twist 20 degrees left, and side bend 15 
degrees right) can be accomplished in one function with three 
loop iterations. We assumed for the model that the maximum 
vertebral joint limit in one dimension will not affect the joint 
limits of another dimension. 

The model includes the spine's rest position because it is the 
position of highest comfort and stability. If the spine is unbend- 
ing in one dimension of movement, It will move towards that 
position of highest comfort in that rotational dimension. What 
determines how much each vertebra bends as the spine moves? 
First, consider one-dimensional rotation, then apply the princi- 
ples to the 3D model. Within one dimension are several factors: 

• Current position 

• Target position 

• Direction of movement: unbending or bending 

• Position of the vertebra 

• Initiator joint 
« Resistor joint 

• Amount of resistance 

• Whether this joint is frozen or participating 

• Whether the position is calculated past the joint limit If so, 
set the position to the joint limit 

Participation of the spine 

A participation vector is derived from the spine's current 
position, target position, and maximum position. This global 
participation represents a 3D vector of the ratio of spine move- 
ment to the maximum range of movement Participation is used 
to calculate the joint weights. 

The following formulas are defined in each of three degrees 
of freedom: 



let 

Target * spine target position . 
Current - spine current position 
Max a spine sum of joint limits 
Ke$(= spine sum of joint rest positions 
participation 

then if spine is bending 

p T foiggr ~ Current 
Max -Current 

else if spine is unbending 

^ „ Tar get - Curren t 
Rest -Current 

Ma/cb 19SU 



Calculation of Joint weights 

The joint positions of the entire spine must sum up to the 
target position. To determine bow much the joint participates,, 
a set of weights is calculated for each joint. The participation 
weight is a function of the joint number, the initiator joint, and 
the global participation derived above. Also, a resistance 
weight is based on the resistor joint, degree of resistance, and 
global participation. To calculate the weight for a given joint, 

let 

i= joint number 

/= joint position 

limit = the joint limit 

rest = the joint's rest position 

w -weight 

p = participation weight 
r = resistance weight 



then if spine is bending 

wpt-rt- (limiti-Jd 

else if spine is unbending 

Wf=Prn*{resu-jb 

The range of weights is from 0 to 1 . A weight of 1 specifies that 
the movement will go 100 percent of the differential between 
the current position and either the joint limit (for bending) or 
the joint rest position (for unbending). A weight ofO means that 
thejoint will moveOpercent of the differential (none at all). The 
weights are a function of the input parameters and global par- 
ticipation. 

Resistance 

To understand resistance, divide the spine into two regions 
split at the resistor joint The region of higher activity contains 
the initiator. Label these regions active and resistive. The effect 
Of resistance is that joints in the resistive region will resist par- 
ticipating in the movement specified by the parameter degree 
pf resistance. Also, joints between the initiator and resistor will 
have less activity, depending on the degree of resistance. 

Resistance does not freeze any of the joints. Even at 100 
percent resistance, the active region will move until all joints 
reach their joint limits. Then, if there is no other way to satisfy 
the target position, the resistive region begins to participate. 

Calculation of joint positions 

If the desired movement is from the current position to one of 
two maximally bent positions, then the weights calculated 
should be 1.0 for each joint participating. The algorithm inter- 
polates correctly to either maximally bent position It also Inter- 
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Figure x 15-14. Wireframe polygonal body pcrforow circular motion. 



pola I us correctly to I he posf lion of highest comfort. To calculate 
the position of each joiol after movement succeeds. 

let 

/a joint position 
/ % new joint position 
Target * spine target position 
Current * spine current position 

M » Target - Current - incremental movement of the spine 



then 
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Bendspine function 

Given the inputs 

• initiatorjoint, 

• rotational movement vector, 

• resistor joint, and 
» resistance 

and a database of joint positions and limits, I he Bendspine func- 
tion outputs the new position of each vertebra after movement. 

Movespine function 

This function drives spinal movement. It captures the joint 
limits and positions nnd» after prompting for input parameters, 
calls Bendspine. Upon returning the new joint positions, Jack 
sets the joint angles and redraws the window. 

Design for a nfew torso 

First, we created a torso using 17 segments corresponding 
with 1 8 vertebral join is. Hacli segment Is (arhitrurily) a hexugo- 
nat slice, with sites (points of attachment) located at the poste- 
rior side of the torso. These sites correspond to the spinal 
vertebral joints. This torso serves as an approximation of 
human anatomy to display movement efficiently. 

Next, we converted more detailed human models based on 
bioslcrcomctrically scanned real (and hence anatomically cor- 
rect) subjects to bodies with 17 contoured slices replacing the 
origlnnlone-scgmenl torso. (Original subject data was supplied 
by Kathleen Robinette of the Armstrong Aerospace Medical 
Research Laboratory at Wright -Patterson Air Force Base.) 

IEP.H Computer Graphics & Applications 
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I'lgure 15. Krc&ihirtg torso Indicates Movement over a 5,0-sce- 
ond hrcalh cycle. Green is 0.0 seconds, red fc J .6 seconds, yellow 
is 4.0 second y. 



These slices were designed to overlap, preventing gaps £rom 
showing up as the torso bends. Each slice was ihcn tiled and 
tested in the Movespine and Bendspinc functions. The final 
result is a realistic torso (hat moves according to the kinematic 
properties ofihc spine. 

Examples 

Our examples show n range of postures and motioit sequences 
involving a hendahle torso and some. The acid test of the/model Is 
real -lime and video animation. We are conlitiulngpur effortslo- 
wufd realistic huhVanpeiformance. 'W " * .* . 

Fign res 1 through 6 show hui nap models based on the! 95th-pe/- 
centile-heigh tmale of a sample of btqsiereometJicaUy scanned real 
subjects. The torsos are displayed In the Jade software system and 
bent using the interactive functions Movcspine and Bcndsplne. 
'J be target state (v,j\ z) indicates the sum of joint positions in three 
degrees of freedom for Ike lS-joitit spine. 

In l-igurcs 7 through 1 2, after t he zero interpolation flag is set 
to "no," the spine is rolled forward and backward around its 
local x axis. Note the figure-eight path drown by the site at the 
top of the head, Th ree other vertebral joints are traced through 
the spnee while rolling the spine. 

In Figures I3and 14, while Jack traces the movement of thTec 
vertebral joints, the spine is circled by varying all three degrees 
of freedom hi coordination. The knees and ankJes bend as the 
spine moves backwards. 

March m\ 



In Figure \$ %: a$ the body breathes, an extension of the upper 
torso occurs with flexion of the lower torso. The pattern re- 
verses on exhalation. On a 5.0-sccond breath cycle, the green 
but line Indicates positions atO.Oscconds, red at 1 ,6 seconds, and 
yellow at 4.0 seconds. 

Future work 

' We plan to use the Movespinc function to generate a vocabu- 
lary of torso gestures, postures, and choreography. Several new 
torsos will be designed and Incorporated into Jack. Body link- 
ages will include n flexible neck with seven ccrvkal vertebrae. 
The model and software can discriminate each cervical joint 
without coding changes. Jn fact, we can input any number of 
vertebral joints to the system. Motion involving interactive 
reach 15 will be redefined to in elude the spine model. Joint limits 
of one rotational degree of freedom should be affected by posi- 
tion within other degrees of freedom. In addition, we should 
enhance the model to indicate the interdependence of lateral 
Inclination and axial rotation doe to the obliquity or the poste- 
rior articular facets. 11 Overall, we must attempt not to oversim- 
plify the model of movement. Experiments so far, however, 
have justified the effort put into the model and appear to have 
added considerably to realistic human figure animation. □ 
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